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Abstract. - We theoretically investigate the static structure factor of an interacting Fermi gas near 
the BEC-BCS crossover at large momenta. Due to short-range two-body interactions, we show that 
the structure factor of unlike spin correlations Sf±(q) falls off as 1/q in a universal scaling region 
with large momentum Tiq and large scattering length. The scaling coefficient is determined by the 
celebrated Tan's contact parameter, which links the short-range behavior of many-body systems 
to their universal thermodynamic properties. By implementing this structure-factor Tan relation 
together with the random-phase approximation and the virial expansion theory in various limiting 
cases, we show how to calculate Sfi(q) at zero and finite temperatures for arbitrary interaction 
strengths, at momentum transfer higher than the Fermi momentum. Our results provide a way 
to experimentally confirm one of the Tan relations and to accurately measure the value of contact 
parameter. 



Strongly interacting fermions are an important feature of many physical systems in 
condensed matter physics, nuclear physics and astrophysics. The recent achievements of 
trapping and cooling ultracold fermionic atoms have made them significant in atomic physics 
as well pQ. With the unprecedented controllability of interactions and geometries [1], these 
atomic systems are prototypes of many important theoretical models. The strength of 
the interactions is governed by a single (i.e., s-wave) two-body scattering length, and can 
be tuned by means of Fcshbach resonances, allowing for a systematic exploration of the 
crossover from a weakly coupling BCS supcrfluid to a Bosc-Einstcin condensate (BEC) of 
tightly-bound molecules. Ultracold atoms therefore provide a new platform for investigating 
the intriguing many-body properties of fermions. A rich variety of theoretical predictions 
has emerged, many of which await verification. Theoretical challenges arise especially when 
the scattering length a is comparable to or larger than the inter-particle distance / [2] • In this 
strongly interacting regime, the system has universal scaling properties that depend on I only 
[3j|4]. These universal properties, however, cannot fully be understood using conventional 
perturbation methods, which are reliable only for a <C I [3[6]. Ab-initio quantum Monte 
Carlo (QMC) simulations are very helpful 015], but suffer from the Fermi sign problem in 
many cases. The challenging nature of strong interactions therefore makes exact results very 
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valuable. 

In this Letter, we discuss theoretically an exact relation for the large-momentum be- 
havior of the spin-antiparallel static structure factor S^i(q), showing that it has a simple 
universal power-law (1/q) tail. This power- law behavior was recently confirmed at Swin- 
burne University, using 6 Li atoms near a Feshbach resonance at nearly zero temperature [5] ■ 
Here, we present the theoretical details. 

The relation we consider, hereafter referred to as the structure-factor Tan relation, be- 
longs to the family of exact relations obtained by Tan in 2005 [TU], which link the short- 
range, large- momentum, and high-frequency asymptotic behavior of many-body systems 
to their thermodynamic properties |10H14j . For instance, the momentum distribution and 
rf spectrum fall off as q~ 4 and w~ 5 / 2 , respectively. All the Tan relations are related to 
each other by a single coefficient I, referred to as the integrated contact density or con- 
tact. The contact measures the probability of two fermions with unlike spins being close 
together [IS]. It also links this short-range behaviour to thermodynamics via the adiabatic 
relation, dE/d(—l/a) = h 2 I / '(Airm) , which gives the change in the total energy E due to 
adiabatic changes in the scattering length. The fundamental importance of the Tan relations 
is due to their wide applicability: to zero or finite temperature, supcrfluid or normal phase, 
homogeneous or trapped, few-body or many-body systems. 

The structure- factor Tan relation for S-\\.(q) follows directly from the short-range be- 
havior of the pair correlation function (r) = J dR. (n-j- (R — r/2) nj, (R + r/2)}, which 
diverges as 
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A Fourier transformation of ni / (r — > 0) then leads to 
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where kp is the Fermi wave- vector and ./V is the total number of atoms. On the right hand 
side of the above equation, we have defined t{q) = [kp/(Aq)\ [1 — 4/{iraq)\. Eq. @ holds 
in a scaling region of sufficiently large q near the unitarity limit (a — > ±oo) so that the 
next-order correction in the bracket (oc l/(aq)) is small compared to the leading term of 1. 
However, the momentum q should always be much smaller than the inverse of the range of 
the interaction potential, so that the pair correlation function has a 1/r 2 behavior. Eqs. (1) 
and (2) were discussed briefly in Ref. [5]. 

The power-law tail of 1/q in the structure- factor Tan relation is more amenable for 
experimental investigation than the q~ 4 or u>~ 5 / 2 tail in the Tan relations for the momentum 
distribution or using rf spectroscopy. The fast decay due to the higher power law in these 
latter two cases imposes very stringent signal-to-noise requirements for studying a given 
range of momenta or frequency. Experimentally, the static structure factor S(q) = S^(q) + 
S^(q) can be readily measured using two-photon Bragg spectroscopy on a balanced two- 
component atomic Fermi gas near a Feshbach resonance |16j . In the large-g limit, to an good 
approximation S^(q) ~ 1. One thus can directly determine the spin-antiparallel structure 
factor S^i(q) = S(q) — 1 and verify a simple 1/q asymptotic behavior, as shown in Ref. [|5]. 

The purpose of this Letter is to quantify the size of the universal scaling region where 
the structure-factor Tan relation @ is valid. For this goal, it is desirable to calculate 
S^i(q) at large momentum (i.e., q > kp) at arbitrary interaction strength and temperature. 
This is a challenging theoretical problem. So far, the static structure factor of strongly 
interacting fermions has only been studied using QMC simulations at zero temperature in 
a homogeneous medium [S]. 

In the present work, we first revisit the zero-temperature problem by extending the scal- 
ing result Eq. © to the weak coupling BEC and BCS limits with \kpa\ <C 1. In the BEC 
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regime, Sf±(q) is mainly caused by ground state molecules. In the opposite BCS regime, 
we use a random-phase linear response theory to obtain first the dynamic and subsequently 
the static structure factor. These two limits are then interpolated smoothly to the unitar- 
ity limit. In this procedure, the zero-temperature contact I is determined from a reliable 
theoretical prediction of the ground state energy along the BEC-BCS crossover, via the adi- 
abatic Tan relation. As an alternative, at finite temperatures close to the Fermi temperature 
(T Tp), we develop a quantum virial expansion theory and calculate the structure factor 
quantitatively at high temperature. In this high temperature regime, the structure-factor 
Tan relation ([5]) can be analytically demonstrated, and the temperature dependence of the 
contact is extracted. The interpolation strategy adopted at zero temperature can also be 
strictly examined. 

Weak coupling limit — We start by considering the weak coupling BEC limit with a small 
and positive scattering length, l/(kpa) — > +00. Here the system is a dilute gas of weakly 
interacting bound molecules, consisting of two atoms with opposite spins, with a relative 
wave function vPg (r) ~ cxp[— r / a\/ {s/2i:ar). Neglecting the correlation among molecules at 
kpr < 1, the pair correlation function is given by n^(r) ~ *^ leading to [S], 

[% (<? > k F )] bec ~ ^ arctan ( y ) ' ( 3 ) 

which is essentially geometry and temperature independent. For qa 3> 1, the first two terms 
in the Taylor expansion of the expression agrees exactly with Eq. ([2]), since Ibec / {Nkp) — 
4ir/(k F a) HD]. 

In the opposite BCS limit of l/(kpa) — > —00, we follow a mean-field picture and treat 
the system as a gas of quasiparticles subjected to a (dynamical) mean-field Hamiltonian 
j!7j . T~L = Hq + (4:Trh 2 a/m)(5n^ + 5n±), where T-Lq is the free particle Hamiltonian and 
dn a = 8n a {v,t) are the space- and time-dependent density fluctuations with respect to 
equilibrium. Here, we consider a normal state since the result differs by an exponen- 
tially small amount from that of the BCS superfluid state. The dynamic structure factor 
So-o-'(q,u>;T) = — lmxaa'{q, w; T)/[%(1 — e ~ huj / k BTy^ ma y j-, c calculated via the density re- 
sponse functions within a (random-phase) linear response theory |17j . which are given by, 

Xaa' = X°aa> + Eo-'V" X° aa » i^T? a/ m)x a >" a , ■ 

The free response function x<x<x' can ^ e wr hten in a simple form in terms of a dimcn- 
sionless function g(q,X;T): Xaa'il^'^) = ~ [ m ^F / (^ 2 )}g(q, A; T)Sa- a >, with q = q/kp, 
A = mu I '(fiqkp), and T = T/Tp. The function g(q,X;T) can be found in standard text- 
books [17] . To leading order in kpa, it is straightforward to show that for a homogeneous 
Fermi gas, the static structure factor Sfi (q) = (2/iV) / dujSfi(q, w; T) is given by 

00 

L T+WJBCS n 2 J 1 - exp[-2gA/T] 

—00 

This weak-coupling equation holds for arbitrary transferred momentum. We may rewrite 
it into the form, [%(<?)] flC s - I B cs/(Nk F )[-4/{irqa)]f(q; f )/(4q/k F ), where at T = we 
have defined f(q; T = 0) = (9/n)q 3 J °° dX Re g(q, A; 0) Im g(q, A; 0) and have confirmed that 
^BCsf(Nkp) ~ 4(kpa) 2 /3 |T^. At large momentum, asymptotically, f(q ^> 1;T = 0) = 
1 + 2/(5q 2 ) + 0(q~ A ). Thus, at the leading order of interaction strengths (~ kpa), the 
structure factor in the BCS limit is in agreement with the Tan relation ([2]). We note, 
however, that the dominant 1/q tail in the strongly interacting regime (i.e., the factor 1 in 
the square bracket in Eq. is lost. The large momentum tail comes from pairing effects 
and therefore cannot be captured by the random-phase approximation for linear response. 

Structure factor from an interpolation strategy. — The excellent agreement between the 
large-g expansion of the weak coupling results, Eqs. and and the Tan relation @ 
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Fig. 1: (Color on-line) Zero temperature (T = 0) spin-antiparallel static structure factor of a homo- 
geneous Fermi gas along the BEC-BCS crossover at q = 3fcf, calculated by using the interpolation 
strategy as described in the text (Eq. [SJ. The red dashed and blue dot-dashed line show the asymp- 
totic behavior in the weak-coupling BEC (Eq. [3} and BCS (Eq. [4} regime, respectively. The inset 
compares our interpolation results (lines) with the QMC simulations (symbols) [8]. The excellent 
agreement at large momentum justifies our interpolation idea. 

suggests that we may use the following extrapolation for the calculation of the structure 
factor at q > kp and at arbitrary interaction strengths: 



In doing so, we remove the severe requirements of q 3> &f and g> 1/a for validating the 
Tan relation @. In Fig. 1, we show the resulting static structure factor Sf± (q = 3kp) 
for a homogeneous Fermi gas along the BEC-BCS crossover. We have determined the 
contact I by applying the adiabatic relation to the ground state energy calculated from a 
pcrturbative Gaussian pair fluctuation theory [B]. In the figure, the weak coupling results 
in Eqs. ^ and (0| are also reported respectively by the dashed and dot-dashed lines. By 
applying a local density approximation, we may also calculate the zero temperature structure 
factor of a Fermi gas in a harmonic trap with frequency ujq and length a; lo = ^Jh/mujQ. 

1/3 

Using the ideal Fermi temperature and vector at the trap center, ksTp = (3iV) hu>o and 
kp = (247V) 1 / 6 /a/ lo , we find that the trapped structure factor differs only slightly with the 
uniform one (see the black line in Fig. 2). 

The accuracy of the interpolation may be justified by comparing our results to QMC 
simulations [5], as shown in the inset for q = 2kp, 3/cf, and 5kp. We find a large deviation 
occurring at q < kp on the BEC side (not shown in the inset), where one has to take into 
account the correlations between different bound molecules. However, this is irrelevant to our 
purpose of examining large momentum correlations. As we shall see, a further justification 
of our interpolation scheme can be obtained at finite temperature, where a quantitative 
calculation of the contact and structure factor is feasible by using a quantum virial expansion. 
In this temperature regime, a random-phase linear response theory predicts for a uniform 
gas, 




(5) 





p-4 



Static structure factor of a strongly correlated Fermi gas at large momenta 



-{x-ty _ p -(x+t)- 



}/t. We find that f(q > l;f) = 



where the function u(x) 
l + f/(q 2 ) + 0(q- 4 )- 

Quantum virial expansion at finite temperatures. - We now calculate the structure 
factor at high temperatures for a trapped Fermi gas, by extending a previous virial expan- 
sion theory jTSHJIJj. At high temperatures, the fugacity z = cxp(u/k B T) 1 is a small 
controllable parameter, even in the strongly interacting limit. We may expand the ther- 
modynamic potential il = — k B T\x\Z 1 where Z = Trexp[— (H — /iJ\f)/k B T], in a series of 
powers of fugacity. By defining partition functions of clusters in the 71-particle subspace, 
Q n = Tr n [exp(-H n /k B T)}, we find that tt/k B T = -[zQ 1 + z 1 (Q 2 - Q\/2) + ■■■}. The 

pair correlation function n^(R =(ri + T2)/2, r = ri — r2j = (n-f- (ri) h\. (r 2 )) can be calcu- 
latcd by taking a functional derivative of the thermodynamic potential with respect to the 
generating elements of a source term J dridr 2 $-|-(ri)$4,(r 2 )7i-t- (ri)ria. (r 2 ), i.e., n^'(ri,r 2 ) = 
S 2 i}(i3f 7 , d^)/(5df5'di)\ i) _ a . Tthc virial expansion form of the thermodynamic potential is 
particularly suitable for taking a derivative of to Q n , and hence for applying the Hcllmann- 
Feynman theorem. 

To leading order, it is readily seen that n^(R,r) = z 2 J2 P ^P(-Ep/k B T) |* P (R,r)| 2 , 
where the summation is over all the two-body states ^p(R, r) of a pair made of atoms with 
opposite spins. In a harmonic trap, the center-of-mass and relative motion parts of the 
wavefunction are separable. By taking the Fourier transform of the relative coordinate (r) 
and performing a spatial integral over the center-of-mass part (R), we finally arrive at 
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where (3 = huJo/k B T, and {Non.} represents the same sum as the first term in the bracket, 
but for a non-interacting system. The relative energies E re i iTl and radial wavefunctions 
iprei,n(r) are known |18j : E re i^ n = (2v n + 3/2)?kJo, where u n is given by \/2Y(— v n )/T(— v n — 
1/2) = a ho /a, and ^ P(S j, n (r) = A n {r/a ho ) exp[-r 2 /(4a 2 ho )]T(-v n )U(-v n , 3/2, r 2 /(2a 2 J), 
where A n is the normalization factor, and T and U are respectively the Gamma and confluent 
hypergeometric functions. 

The accuracy of the virial expansion (|7|) may be examined by applying the same expan- 
sion to the equation of state |21j . which turns out to be quantitatively reliable at T ~ Tp. 
It is straightforward to improve the accuracy by including contributions from Q3, Q4, and 
so on. However, considerable insight can already be obtained from the leading two-body 



contribution shown in Eq. (J7J . At short distance, ip 2 e i ■ 



0) ~ 2irAl[l-2r/a- (2v n 



m(r 2 /a 2 hc 



q, we find from Eq. ([7} that, 
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where B{fi) = Y, n exp[-/3(2i/„ + 3/2)]2tt^ an d 
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is the contact predicted by the virial expansion. We thus explicitly confirm the structure- 
factor Tan relation ^ within second order virial expansion theory. 

In the BEC limit, the two-body bound state VVeZ,n(»") — ^0 (»") = °xp[— r/a]/ (\Z2nar) 
dominates the summation in Eq. Q, with energy E re ^ n ~ — e B = — Ti 2 /ma 2 . The virial 



expansion of the number equation gives [z 2 /(e 
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exp(e B /k B T) ~ N/2 at high 
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Fig. 2: (Color on-line) Temperature dependence of spin-antiparallel static structure factor in a 
harmonic trap. The zero temperature result is obtained by applying a local density approximation 
to the homogeneous structure factor. The inset shows a comparison of the leading virial expansion 
predictions (symbols, Eq. [7} with the interpolation strategy (lines, Eq. [SJ where f(q; T) is aver- 
aged over the trap. This excellent agreement therefore justifies the interpolation strategy at finite 
temperatures. 



temperatures. By inserting these results into Eq. ([7]), we correctly recover the structure 
factor in the weak-coupling BEC limit, Eq. At unitarity, the normalization factor A n is 
analytically known and we find that B(fi) = yj2pK/{e + P - e^) 1/2 ~ ^Jk B T '/(nhujo). Thus, 
at high temperatures T > T F where z ~ (N/2)/(e+^ 2 - e^/ 2 ) 3 ~ (1/6)(T/T F )- 3 , the 
contact is given by IvE^oo/iNkp) — (v / 27t/6)(T/Tf) _5//2 , and therefore rapidly decreases 
with temperature. The structure factor is simply S^i(q) = XvE.oo/(^Nq)[l — 4/(naq)]. 
Finally, in the BCS limit, we may also work out analytically that B{fi) = ^2/tto 2 /(e + @ — 
e -fif/2_ Thig leadg tQ x VE ,BCs/{Nk F ) = [V2^(T/T F y 3/2 /24](fc F a) 2 , agreeing exactly 
with our prediction from random-phase linear response theory at high temperatures. The 
resulting structure factor is given by 



i x Zve,bcs 
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In Fig. 2, we present the numerical results of the leading virial expansion for the static 
structure factor (q = 3fcp) of a trapped Fermi gas. The fugacity z is determined by the 
number equation, N = —dfl/dfi, expanding up to the second order virial coefficient. We 
show also the zero temperature result obtained from the local density approximation. The 
temperature dependence of the structure factor is evident, particularly on the BCS side. 
In the inset, we check explicitly the interpolation strategy with f(q;T). For q > 2Uf, 
the extrapolation results (symbols) are indistinguishable to the numerical virial expansion 
calculations (lines). 

Universal scaling region of our Tan relation. ■ - We are now ready to determine the 
scaling region in which we expect that the scaled structure factor Sf±(q) /t(q) becomes 
momentum independent and converges to the contact of the system. Fig. 3 shows the 
scaled structure factors of a trapped Fermi gas at different transferred momenta and at low 
and high temperatures. On the BCS side, the scaled structure factor is rather insensitive to 
the varying momentum. We observe that the scaling limit can be easily reached at the BEC- 
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Fig. 3: (Color on-line) Scaled spin-antiparallel static structure factor at T — (a) and T = Tf 
(b) for a trapped Fermi gas. Inset in (a): a blow-up of the T — scaled structure factor at the 
BEC-BCS crossover regime. In the universal region with large momenta and strong interactions, 
the scaled spin-antiparallel structure factor converges to the dimensionless contact I/(Nkp). 

BCS crossover regime (|l/fepa| < 0.5), at a relatively small momentum q = 2kp. However, 
to access the scaling limit in a broader region of interaction strengths (i.e., 1/kpa < 2), a 
larger momentum (q — 5kp) is necessary. 

Conclusion. — Based on the structure- factor Tan relation Eq. (2), we present a sys- 
tematic study of the static structure factor of an interacting Fermi gas near the BEC-BCS 
crossover at a transferred momentum q > kp. The scaling region of the structure-factor Tan 
relation is clarified. These predictions can be readily checked in future Bragg experiments 
at crossover, either near the ground state or at finite temperatures. 
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